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Abstract 

The modulation effect in the direct detection of supersymmetric Cold Dark Mat- 
ter ( CDM) particles is investigated. It is shown that, while normally the modulation 
effect in the total event rate is small, < 5%, in some special cases it becomes much 
larger. It also becomes more pronounced in the differential event rate. It may thus 
be exploited to discriminate against background. 
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fL, ! 1. INTRODUCTION 

lO . 

In this work we will study the differential modulation effect of the event rate for 
detecting supersymmetric dark matter, i.e. its variation with respect to the energy 
■ transferred to the nucleus, due to the Earth's motion. 

, There is now ample evidence that most of the matter of the Universe is non 

luminous, i.e. dark Q and is composed of two components. One is the Hot Dark 
Matter (HDM) component consisting of particles which were relativistic at freeze 
out, while the other is the Cold Dark Matter (CDM) composed of particles which 
were non relativistic. There are many arguments supporting the fact that the CDM 
is at least 60% ]|]. There are two interesting cold dark matter candidates: i) Mas- 
sive Compact Halo Objects (MACHO 's) and ii) Exotic Weakly Interacting Massive 
Particles (WIMP's). Since the MACHO 's cannot exceed 40% of the CDM compo- 
nent jl], |3|, there is room for an exotic candidate. The most natural one is associated 
with supersymmetry, i.e. the lightest supersymmetric particle (LSP). 

The most interesting possibility to directly detect the LSP |I| is via the re- 
^ | coiling of a nucleus (A,Z) in the process: 

X + (A,Z) -^ x + (A,Zy (I) 

(x denotes the LSP). In the above process only the elastic channel is open since the 
energy of the LSP is too low to excite the nucleus. In computing the event rate for 
the above process one proceeds with the following steps: 

1) Write down the effective Lagrangian at the elementary particle (quark) level 
in the framework of supersymmetry as described in Refs. |l|, |j . 

2) Go from the quark to the nucleon level using an appropriate quark model for 
the nucleon. Special attention must be paid to the scalar couplings, which dominate 
the coherent part of the cross section and the isoscalar axial current, which strongly 
depend on the assumed quark model ||, |(| . 

3) Compute the relevant nuclear matrix elements [||, |?], ||, |j| |l0| using as reliable 
as possible many body nuclear wave functions. 

4) Calculate the modulation of the event rate due to the Earth's revolution 



around the sun J|, II]. 
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There are many popular targets [12, 13, 14 1 for LSP detection as e.g. 19 F, 23 Na, 
27 Al, 29 Si, 40 Ca, 73 ' 74 Ge, 127 1, 207 Pb, etc. 

In a previous paper ||] we computed the modulation effect h, i.e. the oscillation 
amplitude of the total event rate (see below for its precise definition) , by convolv- 
ing with the LSP velocity distribution the event rate, which, among other things, 
depends upon the relative velocity of the LSP with respect to the Earth. Assuming 
a Maxwell Boltzmann distribution |T| of velocities for LSP, we found that h cannot 
exceed the value of 5% which corresponds to small momentum transfer. The actual 
value of h is quite a bit smaller especially for heavy nuclei and relatively heavy LSP 
(m x > 50GeV). We showed that, in some cases, the quantity h may become neg- 
ative suggesting cancellations between the bins that correspond to small and those 
which correspond to relatively large energy transfers. It is, thus, possible that in 
some energy bins the modulation effect can be larger than the the value of h quoted 
above. 

The event rate depends on many parameters since there exist many contri- 
butions to the above process. The most dominant appears to be the coherent con- 
tribution, which arises out of the scalar coupling originating from Higgs exchange or 
squark exchange if there exists mixing between the L and R squark varieties. It can 
also arise from the time component of the vector current originating from s-quark 
and Z-exchange. The latter is favored from the point of view of the couplings but 
it is suppressed kinematically by factors of (3 2 ~ 10~ 6 owing to the fact that the 
LSP is a Majorana particle. Due to its different dependence on the LSP velocity, 
it yields a higher modulation effect. In addition to the coherent part, especially for 
light targets, when the target spin is non zero, one must include the axial current 
(spin matrix element of the nucleus). 

Our main purpose is to calculate the convoluted differential modulation effect H, 
i.e. the ratio of the part of the differential event rate which depends on the position 
of the Earth divided by that which does not (for its definition see below). If one 
considers each of the above mechanisms separately, H depends only on the LSP 
mass and the size of the nucleus. Knowledge of H may not be adequate, however, 
since one needs to know its value in the energy transfer regime where the event 
rate is the largest and hopefully measurable. One might also need the relative dif- 
ferential event rate, i.e. the ratio of the differential rate to the total rate. If one 
considers the above three mechanisms separately, the relative differential event rate 
is independent of the SUSY parameters or the structure of the nucleon. It depends 
on the nuclear structure only mildly through the form factors. So one can make 
quite accurate predictions which depend only on the nuclear size, the mass of the 
LSP and the low energy cutoff imposed by the detector. 

2. EXPRESSIONS FOR THE RATE 

As we have mentioned in the introduction we only need calculate the fraction of 
the differential rate divided by the total rate which is independent of the parameters 
of supersymmetry. Thus we are not going to elaborate here further on these, but 
refer the reader to the literature J|, ||, 15, 16]. For completeness we only give here 



expressions describing the effective Lagrangian obtained in first order via Higgs 
exchange, s-quark exchange and Z-exchange. We will use a formalism which is 
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familiar from the theory of weak interactions, i.e. 



Leff = -^§{(xi7 A 75Xi)</a + (xiXi)J} (2) 

where 

Jx = N lx (f v + f v r 3 + f° Al 5 + fh^)N (3) 

and 

J = N(f s +f l s rz)N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. Note 
that, due to the Majorana nature of the LSP, Xil^Xi = (identically). The param- 
eters fy, f v , f\, f\, fg, f s depend on the SUSY model employed. In SUSY models 
derived from minimal SUGRA the allowed parameter space is characterized at the 
GUT scale by five parameters, two universal mass parameters, one for the scalars, 
niQ, and one for the gauginos, mi/2, as wen as the parameters tan/3, one of Aq, or 
rr?° le and the sign of [i [17|. Deviations from universality at the GUT scale have also 



been considered and found useful pq ]. We will not elaborate further on this point 
since the above parameters involving universal masses have already been computed 
in some models j|, 19] and effects resulting from deviations from universality will 



be published elsewhere pQ] (see also Arnowitt et al in Ref. JH 

The invariant amplitude in the case of non-relativistic LSP can be cast in the 
form U 

\Mf = EfEl ~ m \ +Pt ' Pf \J \ 2 + \J\ 2 + \J\ 2 

m x 

~ /? 2 |J | 2 + |J| 2 + |J| 2 (5) 

where m x is the LSP mass, \Jq\ and |J| indicate the matrix elements of the time 
and space components of the current J\ of Eq. (0), respectively, and J represents 
the matrix element of the scalar current J of Eq. (Q). Notice that | Jo 1 2 is multiplied 
by 1 (the suppression due to the Majorana nature of LSP mentioned above). It is 
straightforward to show that 

|Jo| 2 = AW)| 2 (fv-f l v^P) (6) 

J 2 = ^|F(q 2 )| 2 (7o_/!±_^) 2 (7) 

|J| 2 = 2^Tll^H [/A°o(q) + /AOi(q)] U)\ 2 (8) 
with -F(q 2 ) the nuclear form factor and 

A A 

Oo(q) = S>C/>" q ' Xj '> Oi(q) = E^0>3(j>' iq - Xj (9) 
j=i i=i 

where cr(j), Xj are the spin, third component of isospin (r^p) = \p)) and 

coordinate of the j-th nucleon and q is the momentum transferred to the nucleus. 
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The differential cross section in the laboratory frame takes the form j|] 

% - T^wh^^-W^ + w + m (10) 

where tun is the proton mass, rj = m x /mNA, £ = pi ■ q > (forward scattering) 
and 

cr = — {G F m N f ~ 0.77 x l(T 38 cm 2 (11) 
The momentum transfer q is given by 

q = got qo = Pr—— (12) 

l + rj 

Some values of qo (forward momentum transfer) for some characteristic values of 
m x and representative nuclear systems (light, medium and heavy) are given in Ref. 
H). It is clear from Eq. ( |i~2| ) that the momentum transfer can be sizable for large 
m x and heavy nuclei (77 small). 

Integrating the differential cross section, Eq. ([l0|), with respect to the azimuthal 
angle we obtain 

o fl A-2Z 2 2 2 (£/3) 2 2t? + 1 ff0 f %A-2Z 2 2 2 
+ Us ~ Is — j — ) F («o£ ) 2~ (1 + 77)2 (-fo ~ ^ — I — ' ( Uo * ^ 

+ (f° A fi m 2 Foo(u f) + 2f A f 1 A n (0)n 1 (0)F 01 (u e) 

+ (/i0 1 (0)) 2 F 11 (no£ 2 )}}2^ (13) 

Where 

^'^) = g— 0^0) PT' p '" =0 ' (14) 

The total cross section a(uo,P), whch has been studied previously (see e.g. |Q, |H), 
can be cast in the form 

^ (—!L\ 2 ^_^^ 5 A 2 \\R 2 ( #0 ! A — 2Z 2 

= (H? {A 4— ) 

+ ifs - /i^^) a ]/o(«o) - y (^(/v - /^) 2 ^W] 

+ (/^ Q (0)) 2 Ioo(%) + 2/i/ifio(0)n 1 (0)J i(«o) 

+ (/AO^O)) 2 /!!^)} (15) 
The quantities I p entering Eq. (^) are defined as 

I p (u ) = (1 + p)t^ (1+p) f ° x 1+ p \F(x)\ 2 dx, p = 0, 1 (16) 

io 

where -F(uo£ 2 ) the nuclear form factor and 

n = q 2 b 2 /2 (17) 
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The integrals I pp >, with p, pi = 0,1, result by following the standard procedure 
of the multipole expansion of the e _iqr in Eq. (||). One finds 

(A,*), 



v(«o) = 2y o edeE— ^ ^o)— • (18) 



For the evaluation of the differential rate, which is the main subject of the present 
work, it will be more convenient to use the variables (v, u) instead of the variables 
(v, £). Thus we get 

+ (/°Oo(0)) 2 F oH + 2f° A f 1 A n (0)n 1 (0)F 01 (u) 

du 
2(/v&) 



+ (fXn^O^iu)}^—^ (19) 



u = g 2 6 2 /2, »r = ^~ (20) 

where p r is the reduced mass and the quantity u is related to the experimentally 
measurable energy transfer Q via the relations 

Q = Qou, Qo = — ^ (21) 

Let us now assume that the LSP is moving with velocity v z with respect to the 
detecting apparatus. Then, the detection rate for a target with mass m is given by 

dN p(0) m , , , 

R = — = \v z \a{u,v) 22 

at m Attin 

where p(0) = OSGeV/cm 3 is the LSP density in our vicinity. This density has to 
be consistent with the LSP velocity distribution (see next section). 
The differential rate can be written as 

dR = P(0) A m \v z \da(u,v) (23) 
m Attin 



where da(u,v) is given by Eq. ( |H 

3. CONVOLUTION OF THE EVENT RATE 

We have seen that the event rate for LSP-nucleus scattering depends on the 
relative LSP-target velocity. In this section we will examine the consequences of 
the Earth's revolution around the sun (the effect of its rotation around its axis will 
be negligible) i.e. the modulation effect. This can be accomplished by convoluting 
the rate with the velocity distribution. Such a consistent choice can be a Maxwell 
distribution |l| 

/(</) = (V^orV^) 2 (24) 
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provided that 

VQ = y/(2/3)(u 2 ) = 220Km/s (25) 

For our purposes it is convenient to express the above distribution in the laboratory 
frame, i.e. 

/(v, v E ) = (V^v y 3 e-( v+v ^ 2 / v2 o (26) 

where ve is the velocity of the Earth with respect to the center of the distribution. 
Choosing a coordinate system in which £2 is the axis of the galaxy, £3 is along the 
sun's direction of motion (vo) and xi = X2 x X3, we find that the position of the 
axis of the ecliptic is determined by the angle 7 ~ 29.80 (galactic latitude) and the 
azimuthal angle uj = 186.3° measured on the galactic plane from the X3 axis (5|. 

Thus, the axis of the ecliptic lies very close to the £2^3 plane and the velocity 
of the Earth is 

v_b = vo + vi = vo + vi( since ~k\ — cosacos~/X2 + cosasin'yx-3 ) (27) 
Furthermore 

COS Ct 

vo • vi = vqv\—= rj vqvi sin"/ cos a (28) 

V 1 + coi 2 7 cos 2 uj 

where vq is the velocity of the sun around the center of the galaxy, v% is the speed of 
the Earth's revolution around the sun, a is the phase of the Earth orbital motion, 
a = 2ir(t — ti)/Tg, where t\ is around second of June and = lyear. 
The mean value of the differential event rate of Eq. (p3|) , is defined by 



/dR\ p(Q) m r I da(u,v) ^ 

\ du I m x AnriN J ' 2 du 

It can be more conveniently expressed as 

(f ) = P -^^-J^)(f) (30) 
\du/ m x Am N v du 

where 

dS f \v z \ da(u,v) 3 

W = i 7py /(v ' vs) ^^ dv (31) 

Thus, taking the polar axis in the direction vg, we get 

A = * fVd« f 1 me -(WB+*»*t)rt dff ("> v ) (32) 
W VfaviJo /-i 1 ^ du y ' 



or 



(-) = -jL- rvdvF Q ( 2 ^)e-^<y^ da ^ V) (33) 
du v 67rv E Jo Vq du 



with 



Fq{x) = Xsinhx - coshx + 1 (34) 
Introducing the parameter 

5 = — = 0.27, (35) 
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expanding in powers of 5 and keeping terms up to linear in it we can write Eq. ( |33| 



as 



+ (/SS!o(0)) 2 Foo(«) + 2/S/Afi„(0)Sl 1 (0)F„ 1 (t.) 

+ (/ini(0)) 2 F„(«)}o 2 (36) 
with 0o = vo/c and 

V2fj, r bvo ^ 



The quantities Fo, Fi, Foo, Fqi, Fli are obtained from the corresponding form factors 
via the equations 

F (n) = F 2 (w)[$E, 0) (aVu) + 0.135 cos a$$\ay/u)] (38) 

F p y(n) = F PiP /(u)[^ 0) (av/u) +0.135 cos aS^ay^)] (39) 

Fi(n) = F 2 (w)[$5 0) (aV«) + 0.135 cos c^VV")] (40) 

= r<iyy 2fe - 1 (exp(- y 2 ))F(2 y )) (41) 
Vo7r Jx 



/6vr 

with Fq(x) given in Eq. (34) and 

r y 2 i 
F 1 ( X ) = 2[(^ + l)coshx-Xsinh X -l\ (42) 

For the cases we considered in this work we find that the quantities F pp i(u) are 
almost the same for all isospin channels. We believe this to be a more general 
result. The value of 0.135 was obtained using sin 7 ~ 0.5 
Combining Eqs. (|(]), (||) and (|38| ) - fl4l|) we obtain 

(*IL\ = B°t°Rr0[l + cosa #(«)] (43) 
\ du> 

In the above expressions R° is the rate obtained in the conventional approach [Bj 
by neglecting the momentum transfer dependence of the differential cross section, 
i.e. by integrating Eq. ([}(]) after the form factors F entering Eq. ( |36|) have been 
neglected. The parameter t is the additional factor needed when the form factors 
are included and the total event rate is convoluted with the velocity distribution. 
RrO is the relative differential rate, i.e. the differential rate divided by the total 
rate, in the absence of modulation, i.e. 

1 dr® 

RrQ=jL^L— (44) 

Note that in the above expressions t° was defined so that the quantity RrO is 
normalized to unity when integrated from u m { n to infinity. From Eqs. ([^) - (^) we 
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see that if we consider each mode separately the differential modulation amplitude 
H takes the form 

H(u) = 0.13 5 |^ (45) 

Thus in this case H depends only on u and a. This means that, if we neglect the 
coherent vector contribution, which, as we have mentioned, is justified, H essentially 
depends only on the momentum transfer, the reduced mass and the size of the 
nucleus. 



Integrating Eq. (43) we get 

R = R°t°[l + cosa h(u , Q min )\ (46) 

where Qmin is the energy transfer cutoff imposed by the detector. The effect of 
folding with LSP velocity on the total rate is taken into account via the quantity 
t . All other SUSY parameters have been absorbed in R . Strictly speaking the 
quantity h also depends on the SUSY parameters. It does not depend on them, 
however, if one considers the scalar, spin etc. modes separately. 

Returning to the differential rate it is sometimes convenient, as we will see later, 
to write it in a slightly different form 

/^\ = R°t°{RrO + cosa Rrl) (47) 
Rrl contains the effect of modulation and is given by 

Rrl = (48) 
t° du y 1 

The meaning of RrO and Rrl will become more transparent if we consider each mode 
separately. Thus for the scalar interaction we get R° — > R° sca i ar and 

dri0) ,2p2A,.W0) 



du 
drW 



= a 2 F 2 (u)&o\aVu~) (49) 

0.135 a 2 F 2 (u)$jp(ai/u) (50) 

For the spin interaction we get a similar expression except that R° — > R® p i n and 
F 2 — ► Fppi. Finally for completeness we will consider the less important vector 
contribution. We get R° — > Ry ector and 

^.^(„),^(^-^^±L«, r( ^] (5!) 

±£ = 0.135^»[*<«(^ - ^gfipjfrPW] (52) 

We see that, if we consider each mode separately, RrO and Rrl are independent of 
all the SUSY parameters except for m x . They depend upon the nuclear physics via 
the relevant form factors. 
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4. RESULTS AND DISCUSSION 



The three basic ingredients of the LSP-nucleus scattering are the input SUSY 
parameters, a quark model for the nucleon and the structure of the nuclei involved. 
Experimentally one is interested in the differential rate. In the present work we 
found it convenient to express it in the manner given by Eq. (|43|), i.e. in terms 
of the parameters R°, t°, RrO and the convolution amplitude H. The parameter 
R° contains all the information regarding the SUSY model. It has been discussed 
previously (see e.g Refs. ||, ||) and it is not the subject of the present work. The 
other parameters will be discussed below. One is also interested in the total rate 
(see Eq. ([llf)) . For this, instead of RrO and H, one needs the convolution parameter 
h. 

The parameter t° expresses the modification of the event rate due to its depen- 
dence on the velocity of the LSP and the folding with the LSP velocity distribution. 
The obtained results, which depend on the LSP mass, the nuclear form factors and 
the detector energy cutoff, Qmin, are presented in Tables la and Ha for four nuclear 
targets of experimental interest. 

The obtained results for h, the modulation of the total event rate, are shown in 
Table lb for Pb and in Table lib for some other nuclei of experimental interest. We 
notice that typically h is quite small, < 5%. Quite surprisingly it can become much 
larger for fairly light LSP and large detector energy cutoff. In other words, in such 
cases as the cutoff energy increases the modulated amplitude decreases less than 
the unmodulated one. There seems, therefore, to be a kind of trade off between the 
total rate and the modulation amplitude . Thus the detector imposed cutoffs may 
yield a bonus of sizable modulation effect, if the event rate is still detectable. 

The quantity which the experiments attempt to measure is the differential rate. 
In the present work we found it convenient to work with the relative differential 
event rate with respect to the energy transfer Q, i.e. the differential rate divided by 
the total rate. Instead of Q we found it convenient to express our results in terms 
of the dimensionless parameter u introduced above (see Eq. (p0|)). The parameter 
u is related to the energy transfer by Q = Qq u with Qq given by Eq. (|2l|). 

We focused our attention on the modulation amplitude which is described either 
by the parameter H (see Eq. (|43|)) or by Rrl (see eq.([48|)). Rrl and H are indepen- 
dent of the SUSY parameters and the structure of the nucleon. Rrl mildly depends 
on the nuclear structure, i.e. it depends on the reduced mass of the system, the 
nuclear form factor and the lower energy cutoff imposed by the detectors. H is even 
independent of the nuclear form factor, but somehow it depends on the size of the 
nucleus. 

Summarizing our results we can say the following : 

1. The nucleus 8 2^ 207 1, i- 

In this case Qq = 40 KeV. We considered both the coherent and the spin contribu- 
tion for 

m x = 30, 50, 80, 100, 125, 250, 500 GeV and Q min = 0, 20, 40 KeV 
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employing the harmonic oscillator form factors of Ref. |22[| . Our results are pre- 
sented in Fig. l(a)-(j). Since the parameter H is independent of Qmin-, it is only- 
shown for Qmin = 0. We see that H rises with u and for the same u it decreases with 
the LSP mass. It can become as large as 25 % for light LSP (see Fig. 1(c)). We 
notice, however, that the event rate drops sharply after u = 0.4, i.e. Q = 16 KeV. 
Thus, the most favored region is around u = 0.2 or Q = 8 KeV (see Fig. 1(b)). We 
also see that H is negative at small u and becomes positive as u increases. Notice, 
however, that the event rate is large at low u (see Fig. 1(a)). Hence we have can- 
cellations in the total modulation amplitude. The analogous results for Q m in = 20 
KeV are shown in Figs. l(d)-(e). The latter results are shifted compared to the 
previous ones by Au = 0.125 but they appear otherwise similar. This is misleading 
since it is the result of the normalization adopted (the area under the curves of 
Figs 1(a) and 1(d) is normalized to unity). Notice that the absolute rates are down 
about a factor of 3 from those at Qmin = 0. We see from Table la that the total 
event rates are very much suppressed for Q m in = 40 KeV. Thus if such cutoffs are 
required by the detector, the process is unobservable. We also present results for 
the spin contribution for the isospin (11) channel in Figs. l(f)-(g) for Q m in = 0. 
Our results for Q m in = 20 KeV compared to those of Qmin = show a similar 
trend as those of Figs. l(d-(e) when compared to those of Figs. l(a)-(b). The 
other isospin channels show behavior similar to the (11) channel || Thus we can 
say in general that the differential rate due to the spin contribution falls quite a bit 
slower compared to the coherent rate as a function of u. We also know that the 
total rate shows a similar trend with respect to uq Furthermore, the quantity 
Rrl is a bit broader, which means that the modulation effect is somewhat favored 
in the spin contribution since a broader energy window around the maximum can 
be selected. For purposes of comparison, we present in Figs. 1(h)- (j) the analogous 
results for Q m in = obtained for the less important coherent vector contribution. 
We see that, in addition to the couplings, the LSP velocity distribution favors the 
vector contribution, but this, of course, is not enough to overcome the suppression 
factor Pq (see Eq. (|36|)) due to the Majorana nature of the LSP. 

2. The nucleus 53 / 127 . 

This nucleus is of great experimental interest [^] due to the advantages of the 
Nal detector. In this case Qq = 60 KeV. We show results for the coherent scalar 
interaction employing the harmonic oscillator form factors of Ref. |2^] for 

m x = 30, 50, 80, 100, 125, 250 GeV and Q min = 0, 45 KeV 
Even though for Q m %n = 45 KeV the total rate is suppressed (see Table Ha), for 
the benefit of the experimentalists we will present the corresponding results for the 
differential rate. We do not show the differential rate for m x = 10 since falls off 
too fast as a function of u. So there is no advantage in going to an energy window. 
Our results are shown in Figs. 2(a)-(c) and Figs. 2(d)-(e) for Q m %n = 0,45 KeV 
respectively. Results for Q m in = KeV are also shown in Figs. 2(g)-(h) in the case 
of the spin contribution for the isospin (11) channel. The other channels show a 
similar behavior. The spin form factors were taken from Ref. @. 
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3. The nucleus u Na 23 . 



This nucleus is a part of the same detector as in the previous one. Here Qq = 630 
KeV. Even though for this light nucleus the spin contribution may be relatively 
more important compared to the coherent one, we only considered in this work the 
coherent contribution in a fashion analogous to the Al case discussed below. The 
parameters t° and h are shown in Tables Ha, lib respectively. In this case the 
detector energy cutoff is 8 — 16 KeV. Our results for the differential rate for zero 
energy cut off are similar to those for Al listed below. For Q m in = 16 KeV they 
are shown in Figs. 3(a)-(b). We see that in all cases the differential rate falls off 
real fast as a function of u. This is not surprising since for such a light system the 
momentum transferred to the nucleus cannot be large. 



4. The nucleus i 3 Al 27 . 



A detector with this nucleus has the advantage of very low energy threshold Qmin = 
0.5 KeV. In this case Qq = 480 KeV. Again only the coherent scalar contribution was 
considered. Both harmonic oscillator and Woods-Saxon form factors were tried. The 
difference between them was small. The results presented were obtained with the 
Woods-Saxon form factors with c = 3.07 and ao = 0.519 fm |p3| . The parameters t° 
and h for various LSP masses and cutoffs are given in Tables Ha and lib respectively. 
In our plots we considered the values of m x = 10, 20, 30, 50 GeV. For larger masses 
the results remain unchanged. For Q m in = 0.5 KeV our results for the differential 
rate are shown in Figs. 4(a)-(c). 

4. CONCLUSIONS 

Detectable rates for the LSP-nucleus scattering for some choices in the allowed 
SUSY parameter space are possible 0. Similar results have been obtained in the 



form of scatter-plots by Arnowitt and Nath[18| and more will appear elsewhere 



[20]. Since, anyway, the event rate is indeed very low, one should try to exploit the 
modulation effect, i.e. the dependence of the event rate on the motion of the Earth. 

In the present work, by convoluting the event rate with the LSP velocity distri- 
bution we were able to obtain the annual modulation effect, both for the coherent 
as well as the spin contribution. We were not concerned with the diurnal modu- 
lation since it is undetectable. This was done both in the total rate as well as in 
the differential rate with respect to the energy transferred to the nucleus. For the 
total rate we found it convenient to write our formalism in terms of three factors 
(see Eq. fliq)). The first one, R°, depends on all the relevant SUSY parameters. It 
represents the total event rate, when the velocity dependence of the cross-section 
and the convolution are neglected. The second, t°, is the modification of the event 
rate due to the velocity dependence of the cross-section and the procedure of folding 
with the LSP velocity. The third is the modulation amplitude h. If one considers 
separately each mode (scalar, spin, vector coherent etc.) t° and h depend only on 
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the LSP mass, the nuclear form factors and Q m %n- The parameter t for various LSP 
masses and a number of nuclear systems as a function of various detector energy 
cutoffs is shown in Tables la, Ha. The total modulation amplitude h is shown in 
Tables lb, lib. We see that it is possible to have a modulation effect which is larger 
than the typical value, h < 5%, but in those cases when the total rate is suppressed, 
e.g. for relatively small LSP mass and large Qmin- So detectors with large cutoffs 
should not be offhand considered to be disadvantaged provided that the total event 
rate is detectable. 

In the case of the differential rate, in addition to the factors t° and h mentioned 
above we needed two more factors (see Eq. (^3|)). The relative differential rate RrO, 
i.e. the differential rate divided by the total rate, and the differential modulation 
amplitude H. If one considers separately each mode, H depends only on the reduced 
mass and the size of the nucleus. The differential modulated rate Rrl depends in 
addition on the nuclear form factors. It is negative at small momentum transfer and 
becomes positive as the momentum transfer increases. As a result, h is always less 
than 5% H and tends to decrease in the case of heavier nuclei. This happens because 
,in the case of Qmin = 0, contributions from different regions of the momentum 
transfer tend to cancel. 

Our main result is that the differential modulation amplitude H can become 
quite large as the momentum transfer increases (see Figs. 1(c), 2(c), and 4(c)). 
Our results are very encouraging. Whether this nice feature, however, can be fully 
exploited by the experimentalists will depend on whether they can exploit the energy 
windows around the maximum of Rrl shown in Figs. 1(b), (e), (g), 2(b), 2(e), 3(b) 
and 4(b). The vector coherent contribution, see Figs. l(h)-(j), shows even better 
features , but unfortunately it cannot be utilized, since the total rate R° associated 
with it is suppressed due to the Majorana nature of the LSP. 

In any event we found many circumstances such that the modulation effect, both 
in the total as well as in the differential event rate, may aid the experimentalists in 
discriminating against background. 
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Figure Captions: 

Fig. 1: The relative differential event rate RrO and the amplitudes for modulation 
Rrl and H vs u for the target 82 -P^ 207 (for the definitions see text). The curves 
shown correspond to LSP masses as follows: 

i) Thick solid line •<=>■ m x = 30 GeV. ii) Solid line •<=>■ m x = 50 GeV. iii) Dotted 
line <=^ m x = 80 GeV. iv) Dashed line <=^ m x = 100 GeV. v) Intermediate 
dashed line m x = 125 GeV. vi) Fine solid line <^=> m x = 250 GeV. vii) Long 
dashed line <^=> m x = 500 GeV. If some curves of the above list seem to have 
been ommitted, it is understood that they fall on top of vi). Note that, due to our 
normalization of RrO, the area under the corresponding curve is unity. 

(a) RrO for the scalar contribution and Q m %n = . 

(b) The amplitude Rrl for the scalar contribution and Q m in = . 

(c) The modulation amplitude H, i.e. the ratio of Rrl divided by RrO for Q m in = 

. 

(d) The same as in (a) for Q m in = 20 KeV. 

(e) The same as in (b) for Q m in = 20 KeV. 

(f) The same as in (a) for the spin contribution in the isospin (11) channel. For 
the other isospin channels the results are similar. 

(g) The same as in (b) for the spin contribution in the (11) channel. 

(h) The same as in (a) for the vector coherent contribution. 

(i) The same as in (b) for the vector coherent contribution, 
(j) The same as in (c) for the vector coherent contribution. 

Fig. 2: The same as in Fig. 1 for the target 53/ 127 . 

(a) RrO for Q min = 0. 

(b) Rrl for Q min = 0. 

(c) H for Q min = 0. 

(d) The same as (a) for Qmin = 45 KeV. 

(e) The same as (b) for Qmin = 45 KeV. The style of the curves is the same as 
in Fig. 1. 

Fig. 3: The same as in Fig. 1 for the target i^Na 23 . The curves shown correspond 
to LSP masses as follows: 

i) Dotted line <^=> m x = 10 GeV. ii) Dashed line <^=> m x = 20 GeV. iii) Long 
dashed line <^=> m x = 30 GeV. iv) Fine solid line •<=>■ m x = 50 GeV. 

For LSP masses heavier than 50 GeV the curves cannot be distinguished from 

iv). 

(a) RrO for Q min = 16 KeV. 

(b) Rrl for Q min = 16 KeV. 

Fig. 4: The same as in Fig. 1 for the target 13 Al 27 and Q m in = 0.5 KeV. (a), (b), 
and (c) refer to RrO, Rrl and H respectively. The style of the curves is the same as 
in Fig. 3. 
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TABLE la. The quantity t° for the target 82-P^ 207 - t° takes into account the velocity de- 
pendence of the event rate and the folding with the LSP velocity distribution. It is computed 
for various LSP masses in the allowed SUSY parameter space. The scalar, the vector coherent 
(k=l) as well as the spin contributions are included. In the latter (11), (01) and (11) indicate 
the possible isospin channels. 











LSP 


mass 


in GeV 






Mode 


Qmin(KeV) 


30 


50 


80 


100 


125 


250 


500 









1.23 


0.728 


0.413 


0.316 


0.246 


0.123 


0.0761 


Scalar 






U.4U4 


U.ool 


n 90.Q 

u.zuy 


U.1D4 


n 1 90. 
u.r/y 


U.UDDo 


U.U40O 






40 





2 x 10" 4 


5 x 10~ 4 


7 x 10" 4 


6 x 10" 4 


5 x 10" 4 


4 x 10~ 4 


Vector 





3.349 


1.735 


0.902 


0.671 


0.509 


0.248 


0.151 


Spin 


(11) 





1.57 


1.298 


0.949 


0.793 


0.661 


0.394 


0.266 


Spin 


(11) 


20 


0.082 


0.512 


0.367 


0.344 


0.312 


0.216 


0.155 


Spin 


(00) 





1.45 


1.13 


0.793 


0.655 


0.542 


0.318 


0.213 


Spin 


(01) 





1.51 


1.21 


0.866 


0.719 


0.597 


0.353 


0.237 



TABLE lb. The same as in Table la for the modulation amplitude h. 











LSP 


mass 


in GeV 






Mode 


Q min (KeV) 


30 


50 


80 


100 


125 


250 


500 


Scalar 




20 
40 


0.0295 
0.1543 
0.2525 


0.0151 
0.0774 
0.1598 


0.0054 
0.0401 
0.0991 


0.0022 
0.0292 
0.0784 


-0.0001 
0.0211 
0.0620 


-0.0005 
0.0070 
0.0314 


-0.0059 
0.0013 
0.0177 


Vector 





0.0543 


0.0621 


0.0571 


0.0560 


0.0553 


0.0545 


0.0543 


Spin 
Spin 


(11) 
(11) 




20 


0.0460 
0.1659 


0.0307 
0.0926 


0.9266 
0.0549 


0.0219 
0.0444 


0.0184 
0.0371 


0.0113 
0.0234 


0.0066 
0.0151 


Spin 


(00) 





0.0421 


0.0349 


0.0238 


0.0195 


0.0163 


0.0100 


0.0056 


Spin 


(01) 





0.0440 


0.0369 


0.0252 


0.0207 


0.0174 


0.0107 


0.0061 
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TABLE Ila. The quantity t° for the experimentally interesting targets 53/ 127 ,n Na 23 and 
13 Al 27 (for definitions see Table la). 











LSP mass 


in GeV 








Target 


r) ( is „\r\ 
QminK^eV ) 


1 n 
1U 


20 


30 


50 


80 


100 


125 


250 




I 







2.16 




1.50 


1.04 


0.689 


0.566 


0.469 0.287 








on 
zu 


n n 
U.U 




0.089 


0.170 


0.162 


0.144 


0.127 0.0855 




Scalar 






n n 
U.U 


0.0014 0.0124 0.0198 0.0201 0.0193 0.0150 











2.13 




1.40 


0.960 


0.651 


0.553 


0.473 0.323 




Spin 


(11) 


on 
zu 


n n 
U. U 




0.075 


0.153 


0.167 


0.164 


0.158 0.137 








45 


0.0 


0.0018 0.0288 0.0483 0.0587 0.0674 0.0781 




Na 







2.33 


2.32 


2.31 


2.30 


2.30 


2.30 


2.30 


2.30 




Scalar 




8 


0.454 


1.19 


1.49 


1.69 


1.69 


1.69 


1.69 


1.69 








16 


0.064 


0.570 


0.907 


1.19 


1.19 


1.19 


1.19 


1.19 




Al 







2.32 


2.31 


2.30 


2.29 


2.29 


2.29 


2.29 


2.29 




Scalar 




0.5 


2.11 


2.22 


2.24 


2.25 


2.25 


2.25 


2.25 


2.25 




TABLE lib. The same as in Table Ila for the 


modulation amplitude h. 














LSP mass in GeV 


Target 


Qmin(KeV) 


10 


20 


30 


50 


80 


100 


125 


250 









0.0508 




.0361 


0.0241 


0.0139 


0.0102 


0..0072 


0..0013 


I 




20 


0.0 




0.1298 


0.0734 


0.0426 


0.0331 


0.0258 


0.0126 


Scalar 




45 


0.0 




0.2194 


0.1294 


0.0740 


0.0588 


0.0474 


0.0267 


I 







0.0501 




0.0344 


0.0241 


0.0180 


0.0166 


0.0157 


0.0149 


Spin 


(11) 


20 


0. 




0.1309 


0.0793 


0.0568 


0.0512 


0.0471 


0.0400 






45 


0.0 




0..2215 


0.1402 


0.1018 


0.0910 


0.0809 


0.0630 


Na 







0.0540 


0.0539 


0.0537 


0.0535 


0.0535 


0.0535 


0.0535 


0.0535 


Scalar 




8 


0.1334 


0906 


0.0793 


0.0715 


0.0715 


0.0715 


0.0715 


0.0715 






16 


0.2039 


0.1237 


0.1030 


0.0911 


0.0911 


0.0911 


0.0911 


0.0911 


Al 







0.0538 


0.0538 


0.0596 


0.0534 


0.0534 


0.0534 


0.0534 


0.0534 


Scalar 




0.5 


0.0598 


0.0563 


0.0553 


0.0545 


0.0545 


0.0545 


0.0545 


0.0545 
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